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ABSTRACT 

Spectroscopic and stereoscopic imaging observations of slow magnetoacoustic wave propa- 
gation within a coronal loop are investigated to determine the decay length scale of the slow 
magnetoacoustic mode in three dimensions and the density profile within the loop system. The 
slow wave is found to have an e- folding decay length scale of 20, OOOlggQgkm with a uniform den- 
sity profile along the loop base. These observations place quantitive constraints on the modelling 
of wave propagation within coronal loops. Theoretical forward modelling suggests that magnetic 
field line divergence is the dominant damping factor and thermal conduction is insufficient, given 
the observed parameters of the coronal loop temperature, density and wave mode period. 



Subject headings: MHD 
oscillations — Waves 

Introduction 



Sun: atmospheric motions — Sun: corona — Sun: oscillations — Stars: 



There has been much theoretical work on the 
propagation of magnetohydrodynamic (MHD) 
waves in solar coronal plasmas, and the advent 
of regular space-based observations of the solar 
corona has subjected theory to empirical test. 
The theory of MHD waves in magneti c cylin - 
ders was described in lEdwin fc Roberts! (|l983f ). 
Since then there have been many papers refin- 
ing the theory (see reviews by .Roberts, .200dt 
Nakariakov fc Verwiclit3l2005l) . 



20081: IVerwichte et al.ll2008[ ). 

In general, thermal conduction is found to be 
the dominant mechanism for dissipation of slow 
modes in the corona. There have been a number 
of observations of propagation in coronal struc- 
tures that have been inter preted as manifesta - 



tions of the slow mode. lOfman et al.l (119971) : 



Deforest fc GurmanI (jlQQSl) : lOfman et al.l (|l999[ ) 



There have been many theoretical studies 
examining the propagation of the slow mag- 
netoacoustic mode in coronal loop structures 
jPorter et al."l994'; Nakariakov ct al. 2000; Tsikl auri fc 
2001; Vc rwichte et al...2001; .Klimchuk et al..,200'3 



Roberts 1 20061 ). Theoretical studies investigating 
the damping of the slow mode have concentrated 
on the effects of thermal conduction, compres- 
sive viscosity, optically thin radiation, gravita- 



describe propagating intensity disturbances in po- 
lar plumes which are found to be consistent with 
the theory of propagating slow modes. There 
have also been many observations of coronal 
loops that are consistent with wave-guiding of 
the slow mode along the loop structure. These 

are^ obse rved as low amplitude intensity and ve- 

■^^S'teSi^)^9)Y,cillations located at the base of quiescent 
coronal loop systems (see review bv iDe Moortdl 
20091 and r eferences within) and a dditionally 



Marsh et al.l (l2008l) : IWang et al.l (l2009t) . Recently, 



Mariska fc MuglachI (|2010[ ) describe long period 



shocks (lOfman et al.l 2000 


: Ofman fc Wand 20021: 


De Moortel fc Hood 2003, 


2004; De Moortel et al. 


2004 ; Mendoza-Briceiio et al. 2004 Erdelvi et al. 



oscillations around 10 minutes observed w ithin 
active regions using Hinode, and Wand ( 2011 ) con- 
ducts a review of standing slow modes observed 
in hot coronal loops. 
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The idea of utilizing observations of wave prop- 
agation to determine proper ties of the co ronal 
plasma was first suggested by lUchidal (|l97nl) and 
[Roberts et al.l |l984l ). termed coronal seismology. 
This idea has since been applied to try and deter- 
mine cor onal properties including dam ping mech- 



2. Observations 



(|2009h . the obser- 



anisms (jOfman fc AschwandenI 12002 ) and dissi- 
pative coeffi c ients ( Nakariakov et al. 19991 ) . In 
Marsh et al.l (|2009l ). observations using the So- 
lar Terrestrial Relations Observatoryt (STEREO) 
spacecraft were used to measure the phase speed 
of a slow mode in the corona directly for the first 
time. This measurement allowed the plasma tem- 
perature to be inferred 'seismologically' via the 
propagation speed of the wave through th e coro- 
nal loop plasma. Marsh fc WalshI ( 2009[ ) tested 
these results using spectroscopic observations of 
the coronal loop system from Hinode; a measure- 
ment of the loop temperature and its spatial pro- 
file, provided an independent confirmation of the 
'seismological' results. 

There is a current, and ongoing, debate re- 
garding the nature of propagating disturbances 
in coronal loop reg ions as being the signature of 
flows or waves ( see De Pontieu fc Mcintosh 2010t 



Verwichte et al.l 120101 ). However, it isn't clear if 
the interpretations in each case are investigating 
the same observational phenomenon. Undoubt- 
edly, both waves and flows exist in coronal struc- 
tures, but the question of which is the cause of 
specific observations of intensity propagations can 
be difficult to answer. The observations presented 
in lMarsh etal . (2009); Marsh & Walsh {200^) are 
interpreted as slow magnetoacoustic waves due to 
the clear periodic nature and the correspondence 
of the propagation speed to the sound speed at the 
temperature measured within the loop. 

In this paper, we determine the density profile 
of the loop system using Hinode observations and 
measure the amplitude attenuation of the wave as 
observed from both STEREO spacecraft. This al- 
lows the true damping scale length of the slow 
mode wave to be directly measured in three di- 
mensions for the first time. This has significant 
theoretical consequences as it places quantitive 
constraints on the modelling of slow mode prop- 
agation within coronal loops and the dissipation 
mechanism of the wave. 



As described in iMarsh et al. 
vations were conducted on 2008 January 10, as 
part of the Joint Observing Program (JOP) 200 - 
'Multi-point, High Cadence EUV Observations of 
the Dynamic Solar Corona'. The STEREO data 
consists of high (30 second) cadence images in the 
17lA bandpa ss from the Extreme -Ultraviolet Im- 
ager fEUVI. IWuelser et al.|[200l . The observa- 
tions begin at 12:00 UT and observe an active re- 
gion loop system which supports slow mode wave 
propagation. The left and right panels of Figure [1] 
show context images of the loop system, along 
with time distance imag es (top) extracted a long 
the indicated paths. See iMarsh et al. I ('2009'), for 
detailed information on the STEREO/EUVI ob- 
servations. 

The same coronal loop system is also observed 
using the E xtreme- Ultraviolet Im aging Spectrom- 
eter (EIS, ICulhane et all l2007l) on the Hinode 
satellite. Figure [U center panel, shows an EIS 
raster of the loop system in Fe XII 195A line in- 
tensity. The spectroscopic raster analyzed here 
uses the 2"slit to build up a 180" x 512"rastered 
image beginning at 18:07 UT. This study contains 
24 spectral windows with a selection of lines de- 
signed to study active regions. 

3. Analysis 

3.1. Preparation of the data 

3.1.1. STEREO/EUVI 

The data are corrected for detector bias, fiat 
field, and photometric calibration, using the SEC- 
CHLPREP routines available within the Solarsoft 
database. The pointing of the EUVI images is 
also co-aligned usi ng the procedure described in 
Marsh et all (|2009[ ). 



3.1.2. Hinode/ EIS 

The EIS data are calibrated using EIS_PREP 
within Solarsoft^ with standard corrections for 
dark current, cosmic rays, hot /warm pixels, dusty 
pixels and an absolute calibration is applied to ob- 
tain the data in units of ergs/cm^ /s/sr/ A 

The pixels along the y axis are binned to in- 
crease the signal to noise ratio within the data, 
resulting in 2" x 2" pixels. The emission line pro- 
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Fig. 1.— Left panel: EUVI 17lA context image of the coronal loop system analyzed in iMarsh et all (|2009l ). 
viewed from STEREO A. The time distance image used to measure the amplitude of the wave propagation 
is extracted along the path indicated by the white lines. To highlight the low frequency propagation, the 
time distance images shown are detrended and convolved in the frequency domain with a narrow Gaussian 
of width 0.14 mHz centred on 1.4 mHz. Center: Fe XII 195A intensity image from Hinode/EIS. The loop 
intensity is extracted from the loop system along the indicated path. The iron emission line ratios are formed 
using the mean intensity along the path to determine the electron density profile as a function of distance. 
Right panel: As left panel, but viewed from STEREO B. Note the spatial coordinates for each image are 
the solar coordinates in the plane of the sky relative to each spacecraft 



files are fitted with multiple Gaussians using the 
Solarsoft routine EIS.AUTO_FIT.GEN, and the 
effects on the line centroids due to the tilt of the 
EIS slit and the orbital variation are corrected. 



3.2. Loop density profile 

The density profile of the loop system support- 
ing the slow mode propagation is measured using 
the Hinode/EIS data. The EIS study contains a 
number of spectral windows that incl ude various 
iron emission lines. lYoung et al detail den- 

sity measurements using Fe XII and Fe XIII lines 
that are observed with EIS. Considering the EIS 
study analyzed here, the Fe XII (186.88/195.12)1 
ratio may be used to determine the density pro- 
file as a function of distance along the loop sys- 
tem, as indicated by the path in the center panel 
of Figure [TJ The density is calculated at a con- 
stant temperature of 0.8 MK, consistent with the 



stant tem perature ot U.g iylK, co nsistent witn tne 
results of (Marsh fc WalshI ( l2009|), using the CHI- 



ANTI v6.0.1 atomic database ( Dere et al.Ml997 



20091 ). To reduce the contamination effect of emis- 
sion from background plasma, a constant back- 
ground is subtracted from the data using the in- 
tensity of pixel [14,22], which is located in a min- 
imal emission region adjacent to the loop system. 



3.2.1. The Fe XII (1 86. 88/1 95. 12) A ratio 

Before determining the Fe XII density ratio, 
the fit to the spectral lines, and possible contribu- 
tions to the line profil es, must be consider ed care- 
fully. As described bv I Young et"al] (|2009l ). Fe XII 
186.88A has a known blend with S XI 186.839A. 
The contribution of this line can be measured 
via branching ratios, using the S XI 188.68A line 
that appears within the Fe XI 188A spectral win- 
dow, as the two sulphur lines have the fixed ra- 
tio of 0.34. In the data analyzed bv lYoung et al 



(|2009l ). the contribution of S XI 186. 839 A to Fe XII 
186. 88A is on the order of 5% and was not cor- 
rected. In the current data, it is not possible to 
fit the S XI 188. 68A line above the noise, indicat- 
ing that the contribution of S XI 186.839A to the 
Fe XII 186. 88A line is even more negligible, and 
thus Fe XII 186.88A is fitted with a single gaus- 
sian. Additionally, the weak line at 186. 98A, in 
the wing of Fe XII 186. 88A, is negligible and so 
neglected. 

The Fe XII 195. 12 A fine includes a blend with 
the Fe XII 1 9 5. 18A transition. As mentioned by 
Young et al. I (l2009l) . the relative strength of this 
line is reduced at lower densities, which is con- 
sistent with the type of quiescent loops observed 
here. The Fe XII 195. 18A hne is not found to 
be significant in this loop system, thus the Fe XII 
195. 12A line is fit with a single Gaussian. 
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3.2.2. Checking the density profile with 

(186. 88/192. 39)A and (186.88/193.51) A 

To confirm the density profile results using the 
Fe XII (186.88/195.12)1 ratio, the 195. 12A line 
may be substituted in the ratio by the density in- 
sensitive Fe XII 192.39A, or 193.5lA lines. These 
lines are available in the EIS study, within the 
Ca XVII 192A window. The (186.88/192.39) A 
and (186. 88/193. 51) A density profiles may then be 
compared to confirm the density profile obtained 
using the 195. 12A line. 

3.3. Fitting the amplitude decay 

A detailed description of the time-distance re- 
duction method and analysis of the wave propa- 
gation is given in Marsh et al. ( 2009f) . The ampli- 
tude of the wave oscillation is fitted at each spatial 
location in time distance images (see Figure [T]) , 
from the A and B spac ecraft, using t he Ba yesian 
fitting code described in lMarsh et all (|2008l) . This 
determines the frequency, phase, and amplitude of 
the wave as a function of distance along the loop 
system, in the plane of the sky as viewed from 
each spacecraft. Beginning with the initial point 
of decay in the amplitude profile, the measured 
amplitude decay is fitted with a function of the 
form 

A(s)=/3e-"-+7, 

where A(s) is the spatial profile of the wave am- 
plitude along the loop system, /3 and 7 are ap- 
propriate constants, a is the decay constant, and 
therefore 1/a is the amplitude decay scale length 
of the wave in the plane of the sky. 

4. Results 

4.1. Loop density profile 

Figure [2] shows the intensity ratio profiles along 
the path within the loop system as indicated in 
the center panel of Figure [TJ These profiles show 
the Fe XII (186.88/195. 12)A, (186.88/193. 51)A, 
and (186.88/192. 39)A ratios. The comparatively 
smaller errors of the (186.88/195. 12)A ratio are 
due to the greater intensity of the 195. 12A line. 
The large error in the (186.88/192.39)A ratio at 
cross-section 6 is due to a badly fitted pixel. 

Figure [3^ shows the CHIANTI density versus 
line ratio curves overplotted with the correspond- 
ing points along the loop path, indicating the 



good overlap between the densities obtained using 
the different ratio diagnostics. Figure Wp shows 
the corresponding density profiles along the loop 
path, demonstrating the good agreement between 
the density profiles obtained using the three iron 
density-sensitive ratios. The three profiles show a 
uniform density profile, consistent within the er- 
rors, along the base of the loop system where the 
wave is observed to propagate. 



Fig. 2. — Intensity ratio profiles, along the 
loop path indicated in the center panel 
of Figure [H denoted by (186.88/195. 12)A 
(square j, (186.88/193.51)A (diamond), and 
(186.88/192. 39)A (star). 



4.2. Three-dimensional amplitude decay 
scale length 

Due to the projection effects of the loop inclina- 
tion and the spacecraft line of sight, the measured 
amplitude profiles of the wave are projected onto 
the plane of the sky as viewed from each space- 
craft. The stereo scopic geometry results from 
Marsh et al. (I2OO9I) are used to rectify the mea- 
sured wave amplitude profiles to the true spatial 
scale, parallel to the wave propagation vector. It is 
assumed that the initial decay in both amplitude 
profiles corresponds to the same spatial location 
within the loop system. To ensure the same re- 
gion of the loop is included within the fit to the 
amplitude decay, the same distance along the loop 
sysyem is fitted in the amplitude profiles on the A 
& B rectified distance scales, shown in Figure ID 

Due to the uncertainty in the inclination of 
the loop, there is an uncertainty in rectified dis- 
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Table 1 

Three-dimensional wave decay scale lengths 



Scale length 1/a (km) 



STEREO A 20, 000^*°™ 

STEREO B 27, 0001^ 



tance scale and a corresponding error on the de- 
cay length scale parameter of the fitted amplitude 
profiles, as plotted in Figure SI Also plotted are 
the max/min fitting errors on the max/min recti- 
fied distance scale uncertainties. The uncertainty 
curves in the plot indicate the max/min ranges 
on the distance scale uncertainty and fitting pa- 
rameters. The maximum and minimum quoted 
scale lengths are the error ranges on the combined 
distance scale and max/min values of the fitted 
parameters. The fit to the amplitude decay gives 
the decay scale length of the wave, of the order 
20,000 km, and associated uncertainties, given in 
Table [T] The fitted length scales of the amplitude 
decay, observed from A & B, agree considering 
these error ranges. 

4.3. MHD Forward modelling 

In this section, we compare our observation- 
ally obtai ned damping length with the theo retical 
model of lOe Moortel fc HoodI (|2003l l2004l ). who 
considered damping by thermal conduction, opti- 
cally thin radiation and area divergence. 



Using the ID MHD model of lDe Moortel k Hood 
(|2003l \2004 . we firstly calculate the value of the 
(dimensionless) thermal ratio d and radiation ratio 
r (Eq uations (5) and (6) in De Moortel fc HoodI 
(120041 ) ■ respectively). Using the observationally 



determined parameters of: loop background tem- 
perature of To = 0.83 MK, background density 
UeQ = 10^'^ cm^"^, oscillation period t = 12 min- 
utes and standard coronal values for all other 
parameters, we find d = 0.0063 and r ~ 0.13. 
Note that these values are respectively a quar- 
ter (for the thermal conduction coefficient) and 
double (for the radiation coefficient) of the values 
used in the De Moortel & Hood studies. Hence, for 
these specific parameters, the damping by thermal 
conduction will be reduced, whereas the effect of 
radiation will be enhanced, compared to the con- 
clusions of iDe Moortel fc HoodI (|2003l . 12004) . 



4-3.1. Thermal conduction 

Let us for the moment consider only thermal 
conduction (setting r — 0). Assuming that all dis- 
turbances are expressed in terms of Fourier com- 
ponents, expi(cLii — kz), the following dispersion 
relation can be obtained 



— iuj'^jdk'^ 



ujk^ + idk'^ = 



(1) 



(see Equation (19) oflDe Moortel fc HoodI (|2003l )) 



As in IDe Moortel fc HoodI l^M i'. we can assume 
a fixed frequency u and solve Equation ([1} numeri- 
cally, where the reciprocal of the imaginary part of 
k corresponds to a damping length. This damping 
length has been plotted as a function of the ther- 
mal ratio d in Figure El (Note that in this figure, 
only the damping length for the slow mode is plot- 
ted, whereas the thermal mode has been omitted 
from this graph.) For the value ol d — 0.0063, we 
find a damping length of the order of 1230 Mm. 
Clearly, this value is far too large to be able to 
explain the observed damping lengths of the order 
of 20 Mm (see Section 4.2). 

4.3.2. Radiation 

Let us now turn to optically thin radiation. 
Here we can di rectly use the exp(— rz/7) d amping 
determined bv lDe Moortel fc HoodI t004\ . where 
r = 0.13. (Note that the variable z in this expres- 
sion is a dimensionless quantity.) Using this value, 
an exponentially damped harmonic oscillation was 
constructed (using the observed values for tem- 
perature, density and the perturbation amplitude 
and period) and subsequently forward-modelled, 
using the response function of EUVI A, to ob- 
tain the synthesised emis sion as would be observed 
by STEREO/EUVI (see lDe Moortel fc Bradsha^ 
^008, for details of the forward modelling proce- 
dure). The predicted intensity (in DN/s/pixel) is 
plotted in Figure [6] as a function of loop length. 
Overplotted (dot-dashed line) is the exponential 
damping, using a damping length of the order of 
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Fig. 3. — a) Chianti curves of emission line ratio 
as a function of density (186.88/195. 12)A (solid), 
(186.88/193.51)1 (dashed), and (186.88/192.39)1 
(dotted), overplotted with the resulting data 
points from the loop path in the center panel 
of Figure [1] b) Corresponding density pro- 
files along the loop path where (186.88/195.12)1 
(solid red), (186.88/193.51)1 (dashed green), and 
(186.88/192.39)1 (dotted blue). 



Fig. 4. — a) The amplitude decay of the 
wave, along the rectified three-dimensional dis- 
tance scale, measured from STEREO A (top) and 
STEREO B (bottom) with the fitted exponential 
decay curve (central solid curve), bounding curves 
indicating the error range due to the uncertainty 
in loop inclination (solid bounding curves), and 
error range with combined fitting errors (dotted 
bounding curves). 



1260 Mm. Similar to thermal conduction, opti- 
cally thin radiation would not be able to account 
for the observed damping lengths. As the damping 
length obtained from thermal conduction and op- 
tically thin radiation are roughly similar, we could 
reasonably expect that the combined effect would 
give a damping length of about half of 1260 Mm, 
i.e. about 600 Mm. However, this is still at least 
an order of magnitude larger than the damping 
lengths estimated from the observations. 

4-3.3. Area divergence 

Finally, we look at the effect of area divergence. 
Assuming a general area divergence of the form 



A{s) — 7r(ro + s tan(6'))^ (with tq the radius at the 
base of the loop, 9 the expansion angle and s the 
distance along the loop), lead s to an amplitude de- 
cay o f the form A(s)~^/^ fsee lDe Moortel &: Hood 
20041 ) ■ Since the precise rotational geometry of 
the loop system is unknown, as an estimate of 
ro and we use the values from the bounding 
loop paths plotted in Fig.[Tl where from STEREO 
A ro — 4.1 Mm and 6 = 5 degrees and from 
STEREO B ro = 5.3 Mm and = 27 degrees. 
The synthesised STEREO/EUVI emission of the 
resulting oscillatory signal ( i.e. a harmonic os- 
cillation with a period of 12 min and damping 
~ 1/ VA) is plotted in Figure [71 To estimate a 
damping length, we have fitted an exponential de- 



6 



0.00 0.01 0.02 0.03 0.04 0.05 
Thermal ratio d 







200 400 600 

oop length (Mm) 



800 



Fig. 5. — The damping length in Mm as a func- 
tion of the thermal conduction ratio d as estimated 
from solving Equation ([T]). The dashed lines indi- 
cate d = 0.0063 and the corresponding the damp- 
ing length. 



Fig. 6.— Synthesised STEREO/EUVI emission 
as a function of length in Mm along the loop. The 
oscillatory signal is damped by optically thin ra- 
diation only. The dot-dashed lines correspond the 
exponential fits to the amplitude decay. 



cay to the first two peaks of the signal, as usually 
only these first two peaks are clearly identified in 
the observations. Using the measured values of 
ro and 9, we obtain a damping length of 105 Mm 
for the STEREO A parameters, and a damping 
length of 65 Mm for STEREO B. Ahhough still 
somewhat larger, this is now of a similar order 
of magnitude as the observed damping length, es- 
pecially keeping in mind that thermal conduction 
and optically thin radiation will cause some addi- 
tional damping. 

4.4. Comparison to dissipation theory 

As referenced in the introduction to this pa- 
per, there has been much work on the theoretical 
modelling of slow mode propagation and dissipa- 
tion within coronal loops investigating effects on 
the wave amplitude such as gravitational stratifi- 
cation and magnetic field divergence, along with 
dissipative mechanisms including thermal conduc- 
tion, compressive viscosity, and optically thin ra- 
diation. 

Ofman fc Wangj (j2002l ) analyzed 35 cases of 



damped Doppler observations in hot >6 MK coro- 
nal loops using the SOHO/SUMER instrument, 
interpreted as sig natures of sta nding slow magne- 
toacoustic modes (jWang et al ][2002, 2003). Using 
a ID MHD model incorporating thermal conduc- 
tion, viscosity and radiation, they found that the 



observed rapid damping of the 5-30 minute period 
oscillations can be explained by thermal conduc- 
tion. 



De Moortel fc HoodI (|2003l ) study a one dimen 



sional model including the effects of thermal con- 
duction and compressive viscosity on 5-minute pe- 
riod waves in 1 MK loops. They find that thermal 
conduction is dominant and that the slow mode 
has a minimum damping length due to ther mal 



conduction alone. iDe Moortel fc HoodI (2004) ex- 



tend this model further including the effects of 
magnetic field divergence and gravitational strat- 
ification, and find that thermal conduction and 
field divergence are the dominant factors. 

investigate analytical 



Klimchuk et al 



(120041) 

and numerical one dimensional models of 5-minute 
waves in ^1 MK loops, and find significant damp- 
ing due to thermal conduction and the observing 
instrument response function. The simulations 
are consistent with damping due to classical ther- 
mal conduction alone without additional damp- 
ing mechanisms, or unusually large conduction 
coefficient. Th ey investigate some pre vious ob- 
servations from iDe Moortel et al.l (|2002l ) with an 
estimated intensity scale length based on their de- 
fined detection length. The advantage of the ob- 
servations that we present here are that they are 
essentially observations of a monochromatic wave. 
Thus allowing a direct measurement of the inten- 




netoacoustic waves in coronal loops of tem- 
peratures around 1 MK and higher, described 
above, find thermal conduction is the domi- 
nant damping mechanism. Thermal conduction 
is a valid explanation of the observed damp- 
ing, even for the observed low oscillation fre- 
quencies, in the type of ~6 MK lo ops observed 
by SUMER (jOfman fc Wand 12002^ : this is also 
the case for 1 MK l oops with oscillation peri- 
ods around 5-niinutes (jPe Moortel fc HoodI 12003 : 
Klimchuk et alj 120041 ) . However, the implications 



Fig. 7.— Synthesised STEREO/EUVI A (solid 
lines) and B (dashed) emission as a function of 
length in Mm along the loop. The oscillatory sig- 
nals are damped by area divergence only, where 
viewed from STEREO A the loop has a basal 
radius tq = 4.1 Mm and loop expansion angle 
— 5 degrees, and viewed from STEREO B 
ro = 5.3 Mm and 6 = 27 degrees. The enveloping 
lines correspond the exponential fits of the ampli- 
tude decay with damping lengths of 105 Mm for 
STEREO A, and 65 Mm for STEREO B. 



sity scale length with the geometrical projection 
effects removed using stereoscopic observations. 

Forward modelling techniques have been ap- 
plied to slow mode propagation (see'De Moortel & 
[2OO8: O wen et al.1 ,2009). De Moor tel fc Bradsha;^ 
pOOSh show that observational parameters such as 
phase, damping rate, and period can be affected 
by plasma emission processes, i.e. the shape of 
the ionisation balance for a particular line can 
lead to a reduction in intensity with an increased 
temperature, depending on the value of the back- 
ground temperature of the plasma and the ampli- 
tude of the temperature perturbations. The loop 
system analyzed here was found to have a uni- 
form temperature prof ile at the base of the loops 
(|Marsh fc Walshll2009^ ■ and, consistent with this, 
the ob served waves also have a constant phase 
speed ( Marsh et al. 20091) . The observed values 
of the loop temperature and density, along with 
the parameters of the observed wave propagation, 
were incorporated into the forward modelling de- 
scribed in Section 1431 

The previous studies investigating slow mag- 



for the class of sub-million degree "cool loops" 
with long oscillation periods (~10-minutes) ob- 
served here has not previously been considered. 
In the case of the temperature and density of the 
loop observed here, we find that thermal conduc- 
tion alone is insufficient to account for the ob- 
served damping of the wave, since the relatively 
cool loop temperature and long oscillation period 
lead to very weak thermal conduction; therefore, 
area divergence has the dominant effect on the 
observed decay of the wave. It should be noted 
that magnetic field divergence may explain the 
observed decay of the wave amplitude, however, 
this is a geometrical effect and not a dissipation 
mechanism of the wave energy. 

5. Conclusions 

The results presented here determine the true, 

Bradshafflree-dimensional, dissipation length scale of a 



slow magnetoacoustic mode in the solar corona. 
The decay scale length at the maximum ampli- 
tude of the wave observed by STEREO A & B 



is 20, OOOIg^^^km and 27, OOOI^^^^^km respectively. 
The agreement between the results from both 
spacecraft, considering the uncertainties, places 
a quantitive constraint on theoretical models of 
slow mode propagation and dissipation mecha- 
nisms within the corona. 

Using the spectroscopic diagnostics of Hin- 
ode/EIS, the density profile of the loop system is 
found to be uniform where the wave propagates. 
This coincide s with a uniform t emper ature profile, 
measured bv iMarsh fc WalshI (|2009l ). The three 
different density diagnostic ratios, under test here, 
give consistent density profiles within the uncer- 
tainties. This suggests that they are are suitable 
diagnostics for these types of loop systems, al- 
though, the (186. 88/195. 12)A yields more precise 
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results due to greater signal to noise. The agree- 
ment also implies that the assumptions regard- 
ing negligible blend contributions to the Fe XII 
195. 12A line, described in Section !??^ are valid. 

The MHD modelling and simulated synthetic 
observations presented here show that thermal 
conduction is insufficient to explain the observed 
short decay length of the waves observed with 
STEREO/EUVI. In contrast to previous studies of 
slow mode damping in coronal loops which can be 
explained by thermal conduction, in the case pre- 
sented here, the divergence of the magnetic field 
is found to be most significant, due to the temper- 
ature, density, and oscillation period within the 
observed loop, damping due to thermal conduc- 
tion is decreased and damping due to radiation is 
enhanced to be of the same order. Field line diver- 
gence alone can account for the observed damping 
to within a factor of three. To determine if the 
combined effects of radiation, thermal conduction, 
and field line divergence can fully account for the 
observed damping length a full numerical study 
incorporating all these effects simultaneously will 
be required in future. 
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